Gene Expression In Colorectal Cancer Research Articles

Transcriptome-wide gene expression outlier analysis pinpoints therapeutic vulnerabilities in colorectal cancer.

Multiple strategies are continuously being explored to expand the drug target repertoire in solid tumors. We devised a novel computational workflow for transcriptome-wide gene expression outlier analysis that allows the systematic identification of both overexpression and underexpression events in cancer cells. Here, it was applied to expression values obtained through RNA sequencing in 226 colorectal cancer (CRC) cell lines that were also characterized by whole-exome sequencing and microarray-based DNA methylation profiling. We found cell models displaying an abnormally high or low expression level for 3533 and 965 genes, respectively. Gene expression abnormalities that have been previously associated with clinically relevant features of CRC cell lines were confirmed. Moreover, by integrating multi-omics data, we identified both genetic and epigenetic alternations underlying outlier expression values. Importantly, our atlas of CRC gene expression outliers can guide the discovery of novel drug targets and biomarkers. As a proof of concept, we found that CRC cell lines lacking expression of the MTAP gene are sensitive to treatment with a PRMT5-MTA inhibitor (MRTX1719). Finally, other tumor types may also benefit from this approach.

Analysis of molecular phenotypes in normal mucosa and colorectal cancer in embryonic anatomical parts of the colon

Background: Differences in the embryonic development of the colonic mucosa determine the physiological embryonic-anatomical asymmetry of its structure and can manifest themselves via different molecular phenotypes (expression profiles) of the colon segments. These molecular characteristics are hypothesized to determine differences in the carcinogenesis mechanisms and influence the prognosis of right- or left-sided colorectal cancer (CRC). Studies of the tumors molecular phenotypes depending on their localization may be of interest for assessment of the prognosis and choice of treatment for CRC.
 Aim: To perform comparative analysis of molecular phenotypes of the normal colonic mucosa and adenocarcinoma CRC tissues depending on the natural embryonic anatomic asymmetry of the colon.
 Materials and methods: We performed a retrospective study of molecular phenotypes (mRNA expression of 61 genes) from different embryonic-anatomical parts of healthy colon and CRC. The normal group included 254 samples of mucosa from three different parts of the colon from 74 healthy donors who had no cancer and no organic abnormalities of the colon, including 90 samples from the right colon, 116 from the left colon, and 48 from the rectum. The CRC group consisted of 154 samples of localized stage T1–4N0–2M0 adenocarcinoma from 154 patients who had not received neoadjuvant radio- and chemotherapy, including 40 samples from the right colon, 54 from the left colon, and 60 from the rectum. The relative mRNA abundance of 61 genes was assessed by reverse-transcriptase polymerase chain reaction. In both groups, the resulting expression phenotypes were compared between the anatomical parts of the colon. Statistical management of the data included the discriminant analysis with stepwise inclusion of variables.
 Results: Based on the assessment of the mRNA level of the studied genes, a discriminant model was built that allows for classification of the normal group samples according to their anatomic origin in the colon with an accuracy of 95.8%. The most significant (p 0.05) for classification are the following 19 genes: CCND1, SCUBE2, TERT, BAG1, NDRG, IL1b, IL2Ra, IL7, ESR1, TGFb, IGF1, MMP9, MMP11, PAPPA, CD45, CD69, TLR2, TLR4, LIFR. The discriminant model built for the CRC group included 27 genes and made it possible to differentiate samples from three parts of colon with an accuracy of 75.2%. A statistically significant (p 0.05) contribution to the samples differentiation by the discriminant model was made by the COX-2, BIRC5, LIFR, TPA, IL1b, MMP11, MMP7, and P16INK4A genes. When combining samples from the two groups into one model in accordance with their embryonic-anatomical origin, there was a clear separation of tumor tissue samples and healthy colonic mucosa in the discriminant function space.
 Conclusion: The analysis of CRC gene expression profiles using the discriminant model showed that genetic changes in the colonic mucosa in CRC flatten the molecular phenotypic boundaries of the embryonic-anatomical parts. These changes are specific to CRC, forming a particular “pathological” molecular phenotype.

Short tandem repeat mutations regulate gene expression in colorectal cancer

Short tandem repeat (STR) mutations are prevalent in colorectal cancer (CRC), especially in tumours with the microsatellite instability (MSI) phenotype. While STR length variations are known to regulate gene expression under physiological conditions, the functional impact of STR mutations in CRC remains unclear. Here, we integrate STR mutation data with clinical information and gene expression data to study the gene regulatory effects of STR mutations in CRC. We confirm that STR mutability in CRC highly depends on the MSI status, repeat unit size, and repeat length. Furthermore, we present a set of 1244 putative expression STRs (eSTRs) for which the STR length is associated with gene expression levels in CRC tumours. The length of 73 eSTRs is associated with expression levels of cancer-related genes, nine of which are CRC-specific genes. We show that linear models describing eSTR-gene expression relationships allow for predictions of gene expression changes in response to eSTR mutations. Moreover, we found an increased mutability of eSTRs in MSI tumours. Our evidence of gene regulatory roles for eSTRs in CRC highlights a mostly overlooked way through which tumours may modulate their phenotypes. Future extensions of these findings could uncover new STR-based targets in the treatment of cancer.

LMTK3 gene expression and the molecular landscape of colorectal cancer (CRC).

172 Background: Lemur tail kinase 3 (LMTK3) plays a critical role in multiple cellular pathways such as Wnt signaling, KIT modulation, and the estrogen receptor pathway. We previously reported that LMTK3 gene polymorphisms are associated with clinical outcome in patients with CRC, and that LMTK3 and estrogen-mediated signaling play a crucial role in CRC tumorigenesis in vitro. Here, we aimed to characterize the molecular features associated with LMTK3 gene expression in CRC. Methods: 20,219 CRC were tested at Caris Life Sciences (Phoenix, AZ) with NextGen Sequencing on DNA (592 genes or WES) and RNA (WTS). Top quartile transcripts per million (TPM) for LMTK3 expression (Q4) were considered high while bottom quartile (Q1) low. Consensus molecular subtypes (CMS) were assessed using RNAseq. Cell infiltration (CI) in the tumor microenvironment (TME) was estimated by QuantiSEQ. X2 and Fisher-Exact tests were used, and significance was determined as p-value adjusted for multiple comparisons ( q < 0.05). Gene expression profiles were analyzed for transcriptional signatures predictive of response to immunotherapy (interferon gamma [IFNγ] score; T-cell inflamed signature [TIS] score). Results: LMTK3 expression was highest in left-sided (0.98 median TPM) followed by right-sided (0.96) and lowest in rectal tumors (0.92). No difference in patient sex distribution was observed between Q1 and Q4 cohorts. Overall, LMTK3 Q1 was associated with TMB-high (≥ 10 Mut/Mb, 11.4 vs 8.7%) and dMMR/MSI-H (8.1 vs 4.6%) ( q < 0.0001); however, these associations were no longer significant in pMMR/MSS tumors. In the pMMR/MSS cohort, LMTK3 expression was highest in CMS4 (1.09 median TPM) followed by CMS2 (0.98), CMS1 (0.88), and lowest in CMS3 (0.70, all intergroup q < 0.05). LMTK3 high showed differential rates of mutations and copy number alterations (CNA) in several genes, including higher mutation rates of APC, TP53, and CDX2 CNA, while lower mutation rates of KMT2A/D, ASXL1, ATM, SMAD4, RNF43 and AMER1 in the overall cohort ( q < 0.05). In the pMMR/MSS cohort, the associations with APC, TP53, SMAD4 and ATM mutations and CDX2 CNA still held true ( q < 0.05). High LMTK3 was associated with higher immune CI in the TME, including B cells, neutrophils, NK cells, CD4+ T cells, monocytes and dendritic cells, while macrophages and Tregs were lower regardless of tumor MSI status (fold change: 0.83-2.4, all q < 0.0001). The IFNγ score was significantly lower in LMTK3 Q4 while the TIS score was significantly higher in LMTK3 Q4 ( q < 0.0001). Conclusions: Our data show a strong association between LMTK3 gene expression and distinct molecular features and TME immune cell infiltration in CRC. These findings suggest that LMTK3 may be an important molecular factor that plays a role in determining the composition of the TME, thus targeting LMTK3 could represent a novel strategy in selected CRC subgroups.

CCR5 and CCL5 gene expression in colorectal cancer: comprehensive profiling and clinical value

BackgroundThe C-C motif chemokine receptor 5 (CCR5)/C-C motif chemokine ligand 5 (CCL5) axis plays a major role in colorectal cancer (CRC). We aimed to characterize the molecular features associated with...

LncRNA VLDLR-AS1 Gene Expression in Colorectal Cancer in Patients from East Azerbaijan Province, Iran

Background: Colorectal cancer (CRC) is the third most common cancer that frequently spreads to other parts of the body, with a low chance of recovery and a high mortality rate. Long non-coding RNAs are considered significant prognostic and diagnostic indicators because they are dysregulated in various cancers and have distinctive expression patterns and high tissue- and cell-specificity. VLDLR-AS1 lncRNA deregulation has been associated with several malignancies. Objectives: This study aimed to evaluate the expression levels of VLDLR-AS1 in CRC. It was the first time this assessment was conducted. Methods: We studied 188 samples, including 94 tumor samples and 94 paired adjacent non-tumor tissues. Total RNA was extracted, and its quantity and purity were assessed. TaKaRa PrimeScript 1st Strand cDNA Synthesis Kit (Kusatsu, Japan) was used for the reaction. The StepOnePlus Real-Time PCR System (Applied Biosystems) was set up to assess the relative expression of VLDLR-AS1. Results: According to the qRT-PCR data, the VLDLR-AS1 expression levels in CRC tissues were significantly lower than in tumor margins (P < 0.0001). In addition, receiver operating characteristic curve analysis demonstrated that VLDLR-AS1 expression can discriminate between tumor and non-tumor samples with the sensitivity and specificity of 72.34% and 51.06%, respectively (P = 0.03, AUC = 0.6274). However, no significant association was found between the expression levels of VLDLR-AS1 and clinicopathological features in CRC patients. Conclusions: The results of this study indicated that VLDLR-AS1 lncRNA is significantly downregulated in the tumor tissues of CRC patients compared to healthy tumor margin tissues. This evidence shows the potential of this gene as a promising biomarker for the early detection of CRC development.

Cooperation of MLL1 and Jun in controlling H3K4me3 on enhancers in colorectal cancer

BackgroundEnhancer dysregulation is one of the important features for cancer cells. Enhancers enriched with H3K4me3 have been implicated to play important roles in cancer. However, their detailed features and regulatory mechanisms have not been well characterized.ResultsHere, we profile the landscape of H3K4me3-enriched enhancers (m3Es) in 43 pairs of colorectal cancer (CRC) samples. M3Es are widely distributed in CRC and averagely possess around 10% of total active enhancers. We identify 1322 gain variant m3Es and 367 lost variant m3Es in CRC. The target genes of the gain m3Es are enriched in immune response pathways. We experimentally prove that repression of CBX8 and RPS6KA5 m3Es inhibits target gene expression in CRC. Furthermore, we find histone methyltransferase MLL1 is responsible for depositing H3K4me3 on the identified Vm3Es. We demonstrate that the transcription factor AP1/JUN interacts with MLL1 and regulates m3E activity. Application of a small chemical inhibitor for MLL1 activity, OICR-9429, represses target gene expression of the identified Vm3Es, enhances anti-tumor immunity and inhibits CRC growth in an animal model.ConclusionsTaken together, our study illustrates the genome-wide landscape and the regulatory mechanisms of m3Es in CRC, and reveals potential novel strategies for cancer treatment.

Value of Coxsackie and Adenovirus Receptor (CAR) Gene Expression in Colorectal Cancer

Value of Coxsackie and Adenovirus Receptor (CAR) Gene Expression in Colorectal Cancer

Ceramide synthase CERS4 gene downregulation is associated with KRAS mutation in colorectal cancer

Ceramide, the central molecule in sphingolipid synthesis, is a bioactive lipid that serves as a regulatory molecule in the anti-inflammatory responses, apoptosis, programmed necrosis, autophagy, and cell motility of cancer cells. In particular, the authors have reported differences in sphingolipid content in colorectal cancer tissues. The associations among genetic mutations, clinicopathological factors, and sphingolipid metabolism in colorectal cancer (CRC) have not been investigated. The objective of this study is to investigate the association between genes associated with sphingolipid metabolism, genetic variations in colorectal cancer (CRC), and clinicopathological factors in CRC patients. We enrolled 82 consecutive patients with stage I–IV CRC who underwent tumor resection at a single institution in 2019–2021. We measured the expression levels of genes related to sphingolipid metabolism and examined the relationships between CRC gene mutations and the clinicopathological data of each individual patient. The relationship between CRC gene mutations and expression levels of ceramide synthase (CERS), N-acylsphingosine amidohydrolase (ASAH), and alkaline ceramidase (ACER) genes involved in sphingolipid metabolism was examined CRES4 expression was significantly lower in the CRC KRAS gene mutation group (p = 0.004); vascular invasion was more common in colorectal cancer patients with high CERS4 expression (p = 0.0057). By examining the correlation between sphingolipid gene expression and clinical factors, we were able to identify cancer types in which sphingolipid metabolism is particularly relevant. CERS4 expression was significantly reduced in KRAS mutant CRC. Moreover, CRC with decreased CERS4 showed significantly more frequent venous invasion.

Inhibition of TRAP1 Accelerates the DNA Damage Response, Activation of the Heat Shock Response and Metabolic Reprogramming in Colon Cancer Cells.

Colorectal cancer (CRC) is one of the major causes of cancer-related mortality worldwide. The tumor microenvironment plays a significant role in CRC development, progression and metastasis. Oxidative stress in the colon is a major etiological factor impacting tumor progression. Tumor necrosis factor receptor-associated protein 1 (TRAP1) is a mitochondrial member of the heat shock protein 90 (HSP90) family that is involved in modulating apoptosis in colon cancer cells under oxidative stress. We undertook this study to provide mechanistic insight into the role of TRAP1 under oxidative stress in colon cells. We first assessed the The Cancer Genome Atlas (TCGA) CRC gene expression dataset to evaluate the expression of TRAP1 and its association with oxidative stress and disease progression. We then treated colon HCT116 cells with hydrogen peroxide to induce oxidative stress and with the TRAP1 inhibitor gamitrinib-triphenylphosphonium (GTPP) to inhibit TRAP1. We examined the cellular proteomic landscape using liquid chromatography tandem mass spectrometry (LC-MS/MS) in this context compared to controls. We further examined the impact of treatment on DNA damage and cell survival. TRAP1 expression under oxidative stress is associated with the disease outcomes of colorectal cancer. TRAP1 inhibition under oxidative stress induced metabolic reprogramming and heat shock factor 1 (HSF1)-dependent transactivation. In addition, we also observed enhanced induction of DNA damage and cell death in the cells under oxidative stress and TRAP1 inhibition in comparison to single treatments and the nontreatment control. These findings provide new insights into TRAP1-driven metabolic reprogramming in response to oxidative stress.

Characterization of MCL-1 in patients with colorectal cancer (CRC): Expression, molecular profiles, and outcomes.

3085 Background: Myeloid cell leukemia 1 (MCL-1) is a member of the BCL-2 protein family and is anti-apoptotic/pro-survival in function. Dysregulation of MCL-1 expression has been reported in several solid tumors, including lung and breast cancer. In CRC, MCL-1 has been associated with resistance to chemotherapeutic drugs and multi-kinase inhibitor regorafenib. Our study aimed to characterize the molecular features associated with MCL-1 gene expression in CRC. Methods: 28,576 CRC samples were analyzed by Caris Life Sciences (Phoenix, AZ) with WTS (Illumina NovaSeq) and NextGen DNA sequencing (NextSeq, 592 Genes and NovaSEQ, WES). MCL-1 expression was stratified by quartiles where top quartile transcripts per million (TPM) were considered high (Q4) and bottom quartile low (Q1). Cell infiltration (CI) in the tumor microenvironment (TME) was estimated by RNA deconvolution analysis using QuantiSEQ. Interferon-gamma and T-cell inflamed signatures were also calculated from RNA data. X2 and Fisher-Exact tests were used, and statistical significance was determined as a P-value adjusted for multiple comparisons ( q < 0.05). Real world survival was obtained from insurance claims data and Kaplan-Meier estimates were calculated for molecularly defined patients. Results: MCL-1 expression was significantly increased in rectal tumors compared to other primary tumor sites (median TPM [mTPM] 20.8 vs 18.9 vs 18.7, in rectal vs left-sided vs right-sided, respectively, q < 0.0001). MCL-1 high tumors were associated with increased rate of PD-L1 IHC levels (Q1 vs Q4: 4.2% vs 9.7% [all CRC]; 3.7% vs 8.9% [pMMR/MSS], q < 0.0001), while MCL-1-low showed marginally increased rates of TMB-high (8.0% vs 7.0% [all CRC]; 3.3% vs 2.6% [pMMR/MSS], P < 0.05). MCL-1-high expression was associated with a higher T-cell inflamed signature and IFN score ( q < 0.0001). MCL-1 high was associated with higher mutation rates of CDKN2A, BRCA2, KRAS and GNAS, while lower mutation rates of TP53, PIK3CA, PTEN, BRAF, APC, FBXW7, AMER1, SOX9, and SMAD2 and copy number amplifications in several genes ( q < 0.05). MCL-1 high TME had higher CI of M1 and M2 macrophages, monocytes, B cells, NK cells, and T-reg infiltration, while dendritic cells and CD4+ T-cell were lower ( q < 0.05). Patients with high MCL-1 expressing CRC had longer survival (Q4 vs Q1: 33.5 vs 24.8 months, HR 0.72, 95% CI 0.67-0.77, P < 0.0001) and longer time-on-treatment with regorafenib although with limited clinical benefit (Q4 vs Q1: 2.03 vs 1.06 months, HR: 0.67, 95% CI: 0.50-0.89, P = 0.005). Conclusions: Our data show a strong correlation between distinct immune biomarkers, TME cell infiltration and MCL-1 expression in CRC. Furthermore, increased tumor MCL-1 expression improved patient prognosis and treatment outcomes. These findings suggest a key clinical role for MCL-1 as an important modulator of anti-tumor immunity and TME and a potential biomarker in CRC.

TCF7L1 Regulates LGR5 Expression in Colorectal Cancer Cells.

Mutations in components of the Wnt/β-catenin signaling pathway drive colorectal cancer (CRC), in part, by deregulating expression of genes controlled by the T-cell factor (TCF) family of transcription factors. TCFs contain a conserved DNA binding domain that mediates association with TCF binding elements (TBEs) within Wnt-responsive DNA elements (WREs). Intestinal stem cell marker, leucine-rich-repeat containing G-protein-coupled receptor 5 (LGR5), is a Wnt target gene that has been implicated in CRC stem cell plasticity. However, the WREs at the LGR5 gene locus and how TCF factors directly regulate LGR5 gene expression in CRC have not been fully defined. Here, we report that TCF family member, TCF7L1, plays a significant role in regulating LGR5 expression in CRC cells. We demonstrate that TCF7L1 binds to a novel promoter-proximal WRE through association with a consensus TBE at the LGR5 locus to repress LGR5 expression. Using CRISPR activation and interference (CRISPRa/i) technologies to direct epigenetic modulation, we demonstrate that this WRE is a critical regulator of LGR5 expression and spheroid formation capacity of CRC cells. Furthermore, we found that restoring LGR5 expression rescues the TCF7L1-mediated reduction in spheroid formation efficiency. These results demonstrate a role for TCF7L1 in repressing LGR5 gene expression to govern the spheroid formation potential of CRC cells.

SETD2 gene expression and the molecular landscape of colorectal cancer (CRC).

184 Background: SETD2, a key methyltransferase modulating histone 3 gene transcription, has been shown to act as a tumor suppressor gene by reducing oxidative stress, colonic inflammation and tumorigenesis in animal models. SETD2 gene expression has been reported to be significantly downregulated in CRC and linked with poorer patient survival. Additionally, SETD2 plays an important role in DNA repair and loss of function mutations have been associated with increased tumor mutational burden (TMB), mismatch repair (MMR) deficiency, and benefit from immunotherapy. Hence, we aimed to characterize the molecular features associated with SETD2 gene expression in CRC. Methods: 15,425 CRC tumors tested at Caris Life Sciences (Phoenix, AZ) with NextGen Sequencing on DNA (592 genes or whole exome sequencing), RNA (whole transcriptome sequencing), and IHC were analyzed. SETD2-high and -low expression were defined as ≥ top and < bottom quartile of SETD2 transcripts per million (TPM), respectively. Cell infiltration in the tumor microenvironment (TME) was estimated by QuantiSEQ. Mann-Whitney U and X2/Fisher-Exact tests were applied where appropriate, with P-values adjusted for multiple comparisons ( q < .05). Gene expression profiles were analyzed for transcriptional signatures predictive of response to immunotherapy (T cell-inflamed) and MAPK pathway activation (MPAS). Real-world overall survival information was obtained from insurance claims data and Kaplan-Meier estimates were calculated for molecularly defined patients. Results: SETD2 expression was significantly increased in left-sided/rectal compared to right-sided tumors (median TPM [mTPM] 55.8 vs 51.1, q < .001). SETD2 mutations were associated with reduced SETD2 expression in pMMR/MSS CRC (mTPM 37.9 vs 54.3 in WT, q < .001), although not statistically significant in dMMR/MSI-H CRC (mTPM 43.6 vs 51.7 in WT, P = .17). Compared to SETD2-high tumors, SETD2-low was associated with increased rates of TMB-high (10.4% vs 8.2%, P = .009), dMMR/MSI-H (7.3% vs 5.6%, P = .02), and PD-L1 IHC levels (4.4% vs 2.5%, q < .002). Only minor differences in gene mutation and copy number rates were observed between SETD2-high and -low tumors , whereas SETD2-high TME was associated with increased immune cell infiltration, including neutrophils, B cells, NK cells, M2 macrophages, and dendritic cells ( q < .05). SETD2-high was also associated with increased T cell-inflamed and MPAS signatures (152% and 159% median increase, respectively, q < .001), regardless of tumor MMR status. Among dMMR/MSI-H CRC, SETD2-high tumors receiving immune checkpoint inhibitors had longer time-on-treatment (HR: 0.39, 95% CI: 0.2-0.76, P = .005) and a trend towards longer OS (HR: 0.29, 95% CI: 0.07-1.2, P = .069). Conclusions: Our data show distinct immune biomarkers and TME cell infiltration associated with SETD2 gene expression in CRC. These findings support a key role for SETD2 in modulating anti-tumor immunity and TME.

Survival of patients with colorectal cancer (CRC) with low expression of homologous recombination proficient (HRP) genes.

225 Background: HR deficient (HRD) CRC has improved outcomes following exposure to DNA damaging agent (DDA) (irinotecan, IR; oxaliplatin, OX) compared to HRP CRC. Low expression of wild type (WT) BRCA1 mRNA is associated with prolonged OS in ovarian cancer; however, this finding has not been investigated in CRC or outside of BRCA. Here, we examine the effect of low expression of HR genes in HRP CRC on post-DDA survival. Methods: 12,860 CRC samples were analyzed by NGS (592, NextSeq; WES, NovaSeq) and WTS (NovaSeq) at Caris Life Sciences (Phoenix, AZ). Samples were classified by RNA expression percentiles. Real world OS was extracted from insurance claims and calculated using Kaplan-Meier estimates for molecularly defined cohorts from first of OX or IR to last contact. Results: Post-IR survival was prolonged with low expression of ATM, CHEK2 and PALB2 in WT ATM, PALB2, and CHEK2 WT CRC, respectively (bottom vs. top 10%, bottom vs. top 25%; p<0.05). Notably, low PALB2 expression (bottom 10%) showed a 14.5-month post-IR benefit compared to high PALB2 expression (top 10%) in WT PALB2 CRC (p=0.003). Post-OX survival was not significantly prolonged with low expression of CHEK2 in WT CHEK2 CRC (bottom 10% vs top 10%, bottom 25% vs top 25%) but was with low expression of ATM (+9.7 months, p=0.02) and PALB2 (+5.6 months, p=0.03) in WT ATM and PALB2 CRC, respectively (all bottom 25% vs top 25%). Conclusions: Here we report the first findings to suggest a novel subclass of CRC defined as the low expression of mRNA of non-mutated HRD genes that exhibit sensitivity to DDA. Low expression of WT ATM, CHEK2, and PALB2 correlates with prolonged OS following IR, post-OX survival was prolonged with low expression of ATM and PALB2. Further characterization defining sensitivity of low expressing HRP genes to DDA may help guide treatment considerations in HRP CRC. [Table: see text]

Lipid Handling Protein Gene Expression in Colorectal Cancer: CD36 and Targeting miRNAs

The reprogramming of lipid metabolism has been highlighted in colorectal cancer (CRC) studies, suggesting a critical role for the scavenger receptor CD36 and fatty acid synthase (FASN) in this malignancy. In this study, we analyzed the gene expression levels of CD36, FASN, the cell surface glypican 4 (GPC4), and the two transporters SLC27A3 and SLC27A4 in 39 paired tumoral and peritumoral tissues from patients with CRC compared with 18 normal colonic mucosae. Moreover, the levels of seven miRNAs targeting CD36 and most of the analyzed genes were evaluated. We found a significant impairment of the expression of all the analyzed genes except GPC4 as well as the differential expression of miR-16-5p, miR-26b-5p, miR-107, miR-195-5p, and miR-27a-3p in the colonic mucosa of CRC patients. Interestingly, CD36 and miR-27a-3p were downregulated and upregulated, respectively, in tumoral tissues compared to peritumoral and control tissues, with a significant negative correlation in the group of patients developing lymph node metastasis. Our results sustain the relationship between CRC and fatty acid metabolism and emphasize the importance of related miRNAs in developing new therapeutic strategies.

Evaluation of Forkhead BOX M1 (FOXM1) gene expression in colorectal cancer.

Forkhead Box M1 (FOXM1)-a key cell cycle regulator is a member of the Forkhead transcription factor family. It plays a key role in embryogenesis and cell proliferation and has been strongly linked to various solid tumors. We sought to understand the regulation of FOXM1 in colorectal cancer (CRC), as well as if and to what extent other clinicopathological characteristics are associated with FOXM1. The investigation comprised 98 CRC samples and normal tissues (controls). All colon cancer patients had a colonoscopy and targeted biopsy. All rectal cancer patients had a CT and MRI. Real-time PCR, Immunohistochemistry, and Western blotting were used to evaluate FOXM1 expression, and the findings were analyzed using SPSS (v.26). FOXM1 mRNA and protein expression were substantially upregulated in tumor tissues, with the majority of these proteins localized in nucleo-cytoplasm. Elevated protein levels of FOXM1 were strongly correlated with lower education level, larger tumor size, lymph node status, lymphovascular invasion (LVI), perineural invasion (PNI), lymph node metastasis (LNM), tumor invasion depth (subserosal and serosal invasion), late stage (III and IV), localization (nucleo-cytoplasmic), intensity (strong) and recurrence. Based on survival analysis, FOXM1 overexpression and nucleo-cytoplasmic localization were associated with shorter disease-free survival while stage and PNI were linked to poorer overall and disease-free survival. According to the results of the Cox regression analysis, stage and PNI were significant predictors of prognosis in CRC patients. FOXM1 expression was elevated in CRC and was linked to reduced disease-free survival. These findings support prior reports and hence FOXM1 can be an important prognostic marker for CRC and a promising therapeutic target. Additionally, we found a link between poor disease-free survival and FOXM1's nucleo-cytoplasmic localization. However, since the sample size of this study was small, further research is needed to validate our findings.

Comprehensive analysis of transient receptor potential channels-related signature for prognosis, tumor immune microenvironment, and treatment response of colorectal cancer.

Transient receptor potential channels (TRPC) play critical regulatory functions in cancer occurrence and progression. However, knowledge on its role in colorectal cancer (CRC) is limited. In addition, neoadjuvant treatment and immune checkpoint inhibitors (ICIs) have increasing roles in CRC management, but not all patients benefit from them. In this study, a TRPC related signature (TRPCRS) was constructed for prognosis, tumor immune microenvironment (TIME), and treatment response of CRC. Data on CRC gene expression and clinical features were retrospectively collected from TCGA and GEO databases. Twenty-eight TRPC regulators (TRPCR) were retrieved using gene set enrichment analysis. Different TRPCR expression patterns were identified using non-negative matrix factorization for consensus clustering, and a TRPCRS was established using LASSO. The potential value of TRPCRS was assessed using functional enrichment analysis, tumor immune analysis, tumor somatic mutation analysis, and response to preoperative chemoradiotherapy or ICIs. Moreover, an external validation was conducted using rectal cancer samples that received preoperative chemoradiotherapy at Fujian Cancer Hospital (FJCH) via qRT-PCR. Among 834 CRC samples in the TCGA and meta-GEO cohorts, two TRPCR expression patterns were identified, which were associated with various immune infiltrations. In addition, 266 intersected genes from 5564 differentially expressed genes (DEGs) between two TRPC subtypes, 4605 DEGs between tumor tissue and adjacent non-tumor tissue (all FDR< 0.05, adjusted P< 0.001), and 1329 prognostic related genes (P< 0.05) were identified to establish the TRPCRS, which was confirmed in the TCGA cohort, two cohorts from GEO, and one qRT-PCR cohort from FJCH. According to the current signature, the high-TRPC score group had higher expressions of PD-1, PD-L1, and CTLA4, lower TIDE score, and improved response to anti-PD-1 treatment with better predictive ability. Compared to the high-TRPC score group, the low-TRPC score group comprised an immunosuppressive phenotype with increased infiltration of neutrophils and activated MAPK signaling pathway, but was more sensitive to preoperative chemoradiotherapy and associated with improved prognosis. The current TRPCRS predicted the prognosis of CRC, evaluated the TIME in CRC, and anticipated the response to immune therapy and neoadjuvant treatment.

Impact of Clinical Data Veracity on Cancer Genomic Research.

Genomic analysis of tumors is transforming our understanding of cancer. However, although a great deal of attention is paid to the accuracy of the cancer genomic data itself, less attention has been paid to the accuracy of the associated clinical information that renders the genomic data useful for research. In this brief communication, we suggest that omissions and errors in clinical annotations have a major impact on the interpretation of cancer genomic data. We describe our discovery of annotation omissions and errors when reviewing an already carefully annotated colorectal cancer gene expression dataset from our laboratory. The potential importance of clinical annotation omissions and errors was then explored using simulation analyses with an independent genomic dataset. We suggest that the completeness and veracity of clinical annotations accompanying cancer genomic data require renewed focus by the oncology research community, when planning new collections and when interpreting existing cancer genomic data.

A bioinformatics framework to identify the biomarkers and potential drugs for the treatment of colorectal cancer.

Colorectal cancer (CRC), a common malignant tumor, is one of the main causes of death in cancer patients in the world. Therefore, it is critical to understand the molecular mechanism of CRC and identify its diagnostic and prognostic biomarkers. The purpose of this study is to reveal the genes involved in the development of CRC and to predict drug candidates that may help treat CRC through bioinformatics analyses. Two independent CRC gene expression datasets including The Cancer Genome Atlas (TCGA) database and GSE104836 were used in this study. Differentially expressed genes (DEGs) were analyzed separately on the two datasets, and intersected for further analyses. 249 drug candidates for CRC were identified according to the intersected DEGs and the Crowd Extracted Expression of Differential Signatures (CREEDS) database. In addition, hub genes were analyzed using Cytoscape according to the DEGs, and survival analysis results showed that one of the hub genes, TIMP1 was related to the prognosis of CRC patients. Thus, we further focused on drugs that could reverse the expression level of TIMP1. Eight potential drugs with documentary evidence and two new drugs that could reverse the expression of TIMP1 were found among the 249 drugs. In conclusion, we successfully identified potential biomarkers for CRC and achieved drug repurposing using bioinformatics methods. Further exploration is needed to understand the molecular mechanisms of these identified genes and drugs/small molecules in the occurrence, development and treatment of CRC.

Genomic Effect of DNA Methylation on Gene Expression in Colorectal Cancer.

Simple SummaryAbnormal DNA methylation is known to regulate gene expression, and its features have been frequently observed in colorectal cancer (CRC) patients. In addition, alterations in DNA methylation can be proposed as biomarkers for cancer prognosis, as they occur in the early stage of carcinogenesis. Although numerous studies have attempted to shed light on the impacts of DNA methylation on gene expression, it is still unclear which specific regions regulate gene expression and how they are associated with patient survival. In this study, we elucidated the intricate relationship between DNA methylation and gene expression. Furthermore, we found genes that were influenced by DNA methylation and were associated with survival; these genes were mainly enriched in immune-related pathways.The aberrant expression of cancer-related genes can lead to colorectal cancer (CRC) carcinogenesis, and DNA methylation is one of the causes of abnormal expression. Although many studies have been conducted to reveal how DNA methylation affects transcription regulation, the ways in which it modulates gene expression and the regions that significantly affect DNA methylation-mediated gene regulation remain unclear. In this study, we investigated how DNA methylation in specific genomic areas can influence gene expression. Several regression models were constructed for gene expression prediction based on DNA methylation. Among these models, ElasticNet, which had the best performance, was chosen for further analysis. DNA methylation near transcription start sites (TSS), especially from 2 kb upstream to 7 kb downstream of TSS, had an essential regulatory role in gene expression. Moreover, methylation-affected and survival-associated genes were compiled and found to be mainly enriched in immune-related pathways. This study investigated genomic regions in which methylation changes can affect gene expression. In addition, this study proposed that aberrantly expressed genes due to DNA methylation can lead to CRC pathogenesis by the immune system.